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Abstract: In previous work we have developed a general method for casting stochastic partial differential equations (SPDEs) 
into a functional integral formalism, and have derived the one-loop effective potential for these systems. In this paper we apply 
the same formalism to a specific field theory of considerable interest, the reaction-diffusion-decay system. When this field 
theory is subject to white noise we can calculate the one-loop effective potential (for arbitrary polynomial reaction kinetics) 
and show that it is one- loop ultraviolet renormalizable in 1, 2, and 3 space dimensions. For specific choices of interaction terms 
the one-loop renormalizability can be extended to higher dimensions. We also show how to include the effects of fluctuations 
in the study of pattern formation away from equilibrium, and conclude that noise affects the stability of the system in a way 
which is calculable. 

PACS: 02.50.Ey; 02.50.-r; 05.40.+j 

Keywords: effective potential, reaction, diffusion, decay 



I. INTRODUCTION 



Geometrical patterns are ubiquitous: From galaxies to living systems, examples abound where a particular spatial 
distribution of some material is preferred versus others out of a seemingly unlimited variety. In many cases, these 
patterns are successfully described by systems of coupled parabolic non-linear partial differential equations. This is the 
case, for example, in chemical kinetics, where such equations summarize the space-time evolution of chemical species 
diffusing and reacting in some confined geometrical region, which makes its presence felt in the boundary conditions for 
the problem. In this way chemical kinetics helps one to understand leopard spots, zebra bands or the radial structure of 
Acetabularia. (See Murray Q, Walgraef ||, and Ball for discussions of many specific examples. Useful background 
references include van Kampen (i) and Gardiner ||.) In these phenomena the values of "reaction constants" play 
a role which reminds one of the role played by coupling constants in determining the vacuum (or ground) state in 
a quantum field theory undergoing spontaneous symmetry breaking, or in the description of phase transitions in 
condensed matter systems. This is not surprising, since in the presence of spatial non-equilibrium patterns [|| "it is 
often the case that the continuous symmetry of the system becomes spontaneously broken" . Because of the above, two 
questions come immediately to mind: (a) what is the effect on an existing pattern of the elimination of fast degrees 
of freedom? and, (b) how do fluctuations affect the stability of an established pattern? These two questions are, of 
course, formally related since both phenomena manifest themselves through noise added to the otherwise deterministic 
equations describing the formation of these patterns. Providing answers to these questions opens the door to the study 
of complex phenomena where either there is no precise explicit knowledge of many of the microscopic details, or else 
where unexpected external perturbations and disturbances show up, and one is nevertheless interested in an explicit 
understanding of the system at long wavelengths. An example is in population ecology, where the presence of illegal 
hunting, or the accidental introduction of some (apparently minor) contaminant, can produce major ecological shifts. 

The most appropriate tool to study patterns of symmetry is, of course, the notion of a potential. If some notion 
of potential is available, then analysis of its extrema leads to the identification of the stable and metastable vacua 
for the system. From there, given the values of the couplings in the system, one determines the ground state for 
the system. This analysis is performed in equilibrium, but the problem in pattern formation, or in reaction-diffusion 
systems, is that these systems are away from equilibrium, often far away from equilibrium, and the usual notion of 
"potential", "vacuum state", and allied concepts are no longer available. Recently, however, we have introduced a 
notion of "effective potential" which generalizes the standard (Quantum Field Theory) notion of an effective potential 
to systems away from equilibrium. The generalization (as discussed below) is such that a very clean and clear parallel 
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can be established with the situation in Quantum Field Theory, and a potential is constructed which has two major 
and distinct pieces: a "classical" contribution and a "fluctuation" contribution. The "classical" contribution plus the 
"fluctuation" contribution determine the "vacuum state" of the system, and therefore this effective potential allows 
the calculation of the effects of the fluctuations on the ground state of the system. The "minimalist formalism" we 
developed in references is to be contrasted with the traditional Martin-Siggia-Rose (MSR) formalism with its extra 
unphysical conjugate fields The minimalist formalism is an extension and outgrowth of the Onsager-Machlup 

approach , and exhibits similarities to the analysis of Crisanti and Marconi |12j . 

The fluctuation dependent piece of the potential involves integrations over the frequency and momentum domains. 
These integrals require the introduction of a cutoff, which, through a Wilsonian-style procedure, leads to a scale- 
dependence of the parameters of the reaction-diffusion-decay system and therefore has an effect on the type of 
instability (and associated patterns) which controls the behavior of the system. The full "effective potential" then 
constitutes a superb tool to incorporate the effects of non-linearities and fluctuations in the patterns produced by 
systems away from equilibrium of the reaction-diffusion type. 

In this paper, after a brief discussion of the "effective potential" in general, and a description of how the noise 
amplitude plays the role of a loop counting parameter, we specialize to the case of reaction-diffusion-decay systems. 
The "potential" is calculated for d = 0, 1,2, and 3 spatial dimensions, and the special case of monomial interactions 
is given separate treatment. We also provide specialized discussions for higher dimensions and general colored noises. 
Then, as an application of the previous results, we study the effects of Gaussian white noise in both Hopf and Turing 
bifurcations by computing appropriate quantities without fluctuations and with fluctuations. We find that the effect 
of noise is to shift the symmetric states of the system, as well as to change the nature of the linear instabilities that 
may occur as perturbations around these new states. We end by offering some conclusions. 

To fix the notation: In a companion paper (?]], we discussed classical field theories subject to additive stochastic 
noise rj(x, t) described by the equation 



D(j>(x, t) = F[<j>{x, t)} + rj(x, t) 



(1) 



Here D is any linear differential operator, involving arbitrary time and space derivatives, that does not explicitly 
involve the field <j>. The function F[(j)] is any forcing term, generally nonlinear in the field <f>. These stochastic partial 
non-linear differential equations (SPDEs) can be studied using a functional integral formalism which makes manifest 
the deep connections with quantum field theories (QFTs). We showed that if the noise is translation-invariant and 
Gaussian, it is possible to split its two-point function into an amplitude A and a shape function g2(x,y), as follows 



Gn(x, y) d = Agi(x 



V) 



(2) 



with the convention that 



J d d x dt g2 l {x,t) = 1 = 



g^(k = 6,iu = 0) 



(3) 



Then the one-loop effective potential for the SPDE is 



(2 7 r) d + 1 
0(A 2 ) . 



In 



(2?t(fc, £ 



<5Ft 



(4) 



Here (f> is any convenient background field. The above result is qualitatively similar to the one-loop effective potential 
for scalar QFT ||-||: 



V[<t>;U = V(</>) 



dfik&UJ 



d+l 



In 



§2V 

8<p 8<f> 



k 2 



(5) 



as can be seen by simply comparing equations (||) and (||). Moreover, as argued in Q, this effective potential for 
SPDEs inherits many of the interesting features of the effective potential for QFTs. In particular, minima of the 
effective potential for a given SPDE correspond to homogeneous and static solutions of the stochastic equations of 
motion (that is, homogeneous and static expectation values of the stochastically driven field 4>, with the averaging 
done with respect to noise realizations) . 
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In a second paper fL7j] , we applied this formalism to the Kardar-Parisi-Zhang (KPZ) equation, obtaining an inter- 
esting ground state structure (including dynamical symmetry breaking). In the current paper we perform a similar 
analysis for the reaction-diffusion-decay system described by the equation: 

^-- l yW 2 U = P -j<l> + P(fa+V- (6) 

This equation can be used, for instance, as a model to describe the dynamics and spatial distribution of the concen- 
tration of a chemical reagent, when it is subject to both diffusion (via v) and decay (via 7). P(4>) is some ultra-local 
function of the concentration (typically a polynomial in the concentration, but not always) and it represents the 
reaction kinetics |l]-|E|. We have included a "tadpole term" Pq since, as we will see, its presence is essential for the 
consistency of the ultraviolet renormalization program. Many examples of reaction-diffusion equations abound in the 
literature. By way of example, we mention just two model equations that the reader may wish to keep in mind. 
One goes under the name of "amplitude equation" , generally complex, resulting from reducing a (non-stochastic) 
reaction-diffusion equation in the vicinity of an instability point: 

d t A = uV 2 A + f(\A\ 2 )A, (7) 

where the diffusion coefficient can be complex: v = vr + ivi. The study of amplitude equations is very useful for 
determining the basic geometry of patterns that can emerge near the instability point |§H. This equation clearly 
involves only one field degree of freedom, but the generalization is straightforward. For several species of, e.g., chemical 
reactant the field (concentration) cj) is promoted to a vector fa(x,t). The diffusion coefficient and decay rate are then 
promoted to matrices, the noise to a vector, and the reaction kinetic function P(fa) to a vector-valued functional with 
tensorial coefficients. 

Si j ^ - ufi V 2 ) fa = (P )i - 7/ fa + Pi(<f>j) + Vi ■ (8) 

Typical examples of this sort involving two components go under the heading of "activator-inhibitor" models of e.g., 
biological pattern formation, a noise-free example of which is provided by the Gierer-Meinhardt mechanism: 

d t A = D A V 2 A + h - k 2 A + , (9) 

B 

d t B = D B V 2 A + k 4 A 2 - k 5 B . (10) 

Here fa — A is the activator and fa = B is the inhibitor. The ki are reaction constants. 

Finally, we note that these systems provide a viable framework for the notion of "self-organizing systems" . 

The remainder of this paper is organized as follows. In Section II, following the general procedure developed in [Q, 
we construct the complete one- loop effective potential for the class of reaction-diffusion equations in (^) . Once this 
is done, we carry out a detailed analysis of the resulting potential as a function of spatial dimension for d — 0, 1, 2, 3. 
Special attention is also paid to monomial interactions, the case of higher spatial dimensions, and the case of correlated 
Gaussian noise. In Section III we turn to a discussion of the impact that noise and fluctuations can have on the onset 
of instabilities and pattern formation. We conclude in Section IV with a discussion of our results. Certain technical 
issues having to do with functional Jacobian determinants are collected in Appendixes A and B. The general Feynman 
rules for the reaction-diffusion equation (|^) are presented in Appendix C. An integral needed in the computation of 
the one-loop effective potential is calculated in Appendix D. Finally, it should be borne in mind that the main body 
of the paper adopts the so-called Stratonovich calculus — Appendix E indicates the changes that are required if one 
wishes to consider the so-called Ito calculus. 

II. EFFECTIVE POTENTIAL: REACTION-DIFFUSION-DECAY SYSTEMS 

We begin this section by computing the one-loop effective potential associated to the class of reaction-diffusion 
equations (||). We first carry out the computations for arbitrary Gaussian noise and arbitrary spatial dimension, 
only later specializing to white Gaussian noise and examining the particular features of the effective potential for 
dimensions d = 0, 1, 2, 3. 

To avoid unnecessary clutter, let us re- write equation (^) as 

^- t ~vy 2 \<j } ^P(fa+r 1 . (11) 
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Here P(<f>) is an arbitrary polynomial in the field 6. Any tadpole contribution Po, as well as any decay term —70, 
have now for convenience been subsumed into P(</>). 

To apply the general analysis provided in [Q, which led to equations (0)-(Q) above, to a homogeneous and static 
background field, 4>{x, t) — constant, we simply make the identifications 



FW\ - P(0); 



SF 

86{x) 



5 2 F 



56{x)56{y) 



P"(<f>) 



(12) 



as follows by comparing equations (|lj) and (11). At tree- level (zero- loop) in the loop-counting parameter A, the 
equations of motion become |?J 



SF' 



F[<j>] = J 



P'{6)P{6) = J 



(13) 



This is a polynomial equation for cf> and therefore has a finite number of roots. In particular, for J = let </>o be one 
of the roots of P'(<f))P(<fi) — 0. This polynomial always has at least one real root. [Proof: let P(</>) be of degree n, 
then P 2 (4>) is of degree 2n, and P'(0)P((/>) is of degree 2n — 1, which is always odd. Thus P' \6)P(6) must cross the 
abscissa at least once, so there must be at least one real zero.] 
The zero-loop effective potential is now 



V, 



zero— loop 



[6; </) O ]-i[P 2 (0)-P 2 (0 o )] 



(14) 



This zero-loop effective potential is both a generalization (because it includes the effects of non-linearities) and a 
specialization (because it treats static fields) of the Onsager-Machlup action for stochastic mechanics jj],[nj[ll) . 

We start the one-loop computation by noting that for the reaction-diffusion-decay system the linear differential 
operator needed is given by 



D-°— = d t - v v 2 ~ P'{6) 



-iw + vk 2 - P'{6) , 



in configuration and Fourier variables. We then have for the adjoint quantity 

<SP f 



so that 




5cj) 



D- 



SF_ 

56 



-d t - ^V 2 - P'{6) 



= -d 2 t + [vV 2 + P'(6)f 



+iuj + vk 2 - P 1 \6) 



W + [vk 2 - P'{6)f 



(15) 



(16) 



(17) 



Using the previous equations, which are specializations of the general ones presented in 
potential is given by 



the one-loop effective 



Vf 



to] = iP 2 (A) + $ A 



d d k duj 

(2TT) d + 1 



In 



g 2 (k^)P{<j>)P"(^) 
uj 2 + [vk 2 - P'{6)] 2 



0(A 2 ) 



Equivalently 



V[S-6 ] = ±P 2 (6) + ±A 



d d kdiv 

W) 



d+1 



hi 



lo 2 + [vk 2 - P'{6)] 2 + g 2 (k, c)P(0)P"(0) 
oj 2 + [vk 2 - P'{S)} 2 



)) + 0(A 2 



(18) 



(19) 



At this stage we should be explicit about some technical details: First, the way we have chosen to treat the functional 
Jacobian is equivalent to choosing the Stratonovich calculus for the stochastic noise |t| . The modifications attendant 
on the choice of the Ito calculus are in some ways a simplification of the current procedure, but in other ways lead 
to additional technical complications. Appendix E discusses some features of the Ito calculus. Second, in deriving 
the formulae above it has been assumed (in performing the functional integrations) that the field 4> is unrestricted, 
and can take on all values from — oo to +oo. Strictly speaking, this is of course not the case if <fi represents a 
concentration, but this is not a serious restriction. One can deal with this either by (1) choosing the forcing term 
P{4>) to strongly suppress negative values of 6, and then taking a suitable limit, or more prosaicly (2) by realizing 
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that in making the one-loop approximation we have already assumed that fluctuations are in some sense small, so 
that if we look at quadratic fluctuations around some positive value of the background field <f>o, then the error made 
in letting S(j) — <\> — 4>q range over the entire real line is a higher order effect [at least 0(A 2 ).] 

This is as far as we can go without making any further assumptions about the additive noise. For instance, one 
standard choice is temporally white, which means delta function correlated in time, so that g2(k,co) — > 92{k) is a 
function of k only. Let us define X 2 = [vk 2 - P' (</>)] 2 + g 2 (k)P(<p)P" (<p) and Y 2 = [vk 2 - P' (</>)] 2 . It is easy to see 
from its definition that Y 2 is real andpositive, and so is Y. If X 2 is positive, we can make use of the standard integral 
identity [X and Y are positive, see flBf , equation (4.222.1)], namely 



dcu In 



Y 2 



= 2n(X-Y) , 



(20) 



to re- write equation (|l9|) as 



V[frM = hP 2 {<P) + hA J ^[vk 2 - P'm 2 + ~g 2 {k)P^)P"^) K 2 - P'(<f>)\\ - (<f> - <f>o) + 0(A 2 ) . (21) 

If X 2 is real and negative, we can define X 2 = — Z 2 , with Z a real positive number. In this case we must make use 
of the following integral (see Appendix D) 



°° &u\A tjj2 ~ Z \t lt \=^{±iZ-Y) , 



2 + Y 2 



(22) 



to get the appropriate result; here e is a real, small, and positive number. 
Averaging over the ± justifies the following prescription 



+ oo 



dcu In 



LU 2 + Y* 



duj Re 



In 



uj 2 + Y 2 



-2irY 



(23) 



This prescription is guaranteed to preserve the reality of the effective potential. In general we should write 



V[^U = ±P 2 ^) + ±A 



d d k 
(2-K) d 
0(A 2 ) , 



Re 



\vk 2 - P'(<P)} 2 + g 2 {k)P{(j>)P"{(t>) 



-\vk 2 -P'(<j>)\ 



(24) 



where we have used the w-integral to recast fll9| ) into (|24]). Remembering that \im^^g2(k) = 1 [recall (||)], it is clear 

that there are no infrared divergences (k — > 0), at least for this effective potential at one- loop order. To investigate 
the ultraviolet (i.e., short-distance) behaviour, it is useful to re-express this in the form 



V[</>;<M = ^ 2 (0) + |-4 



d d k 



vk 2 -P'{4>) 




g 2 (k)P(ct>)P»(4>) 



[vk 2 - P'(4>)Y 



- 1 



{2-kY 

-(cP^^ + OiA 2 ) . (25) 
It is clear now that this effective potential will be finite provided the spatial part of the noise spectrum satisfies 



d d k 



\k 2 + a 



< oo . 



(26) 



Thus, while the noise acts as the source of the fluctuations, it can also serve as the regulator to keep physical 
quantities UV-finite, as should be clear from (p6|). Indeed, whatever modifications one might make to the noise in the 
ultraviolet region will have no consequence for the long- wavelength or hydrodynamic limit. The noise is intended to 
model fluctuations above a certain limiting resolution length/time scale, thereby making its short distance behavior 
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immaterial, in so far as one is interested in studying the long distance, long time asymptotic behavior of the stochastic 
model []. 

For definiteness, let us now take the spatial noise spectrum to be cutoff white, i.e., 

3 2 (£)=£b(|fc|)=e(A-fe). (27) 




With this choice of noise, we can Taylor expand the square root of equation (|2Jj) in the ultraviolet regime to see 
that there is a divergent term proportional to P(4>)P" ((f>) A d ~ 2 , and a subdominant divergent term proportional to 
P{4>) P'(^>) P"{<t>) A d ~ 4 . Since the classical (tree-level) potential is just P 2 {4>), to have any hope of absorbing the 
infinities into the bare action we must have d < 4. That is: the reaction-diffusion-decay system, subject to white 
noise and for any polynomial P(<f>), is one-loop ultraviolet renormalizable only in 0,1,2, and 3 space dimensions. In 
space dimensions the reaction-diffusion-decay system reduces to a Langevin reaction-decay system which is still 
interesting, (see below). In d = 1 the reaction-diffusion-decay system is in one- loop finite. (Strictly speaking the 
claim of one-loop renormalizability also requires the investigation of the wavefunction renormalization; this is beyond 
the scope of the present paper. For related discussion see Gjl.) The assertion that arbitrary polynomial reaction 
kinetics can be renormalizable in low dimensions should not (with hindsight) be alarming. After all, exactly the same 
thing happens for quantum field theories in d = 2 spacetime dimensions, where P(4>)2 is renormalizable for arbitrary 
polynomials. If we restrict the form of the polynomial occurring in the reaction-diffusion-decay system, we can have 
one-loop renormalizable theories in a higher dimensions. We will come back to this particular point later. 
To be more explicit, we expand the unrenormalized one- loop effective potential as follows 

Vfe M = + kA I ^ > : { ( ! ] ; }"(*->*)) + 0(A 2 ) . (28) 

This expansion only makes sense if \P(4>)P" (<j))\ < \vk 2 — P'(0)] 2 for every value of |fc|. This requires both P'(4>) < 
and \P(4>)P" (4>)\ < [P'{4>)] 2 . The relevant integrals, after a rescaling, areg 

InMn, d) ^ r dx ^(1 + x 2 f- 2 - = jB(d/2, 2n - 1 - d/2) = ^ (^2)^ - 1 - d/2) 
Jo ~ 2r(2n-l) 

These integrals converge for n > (rf+2) /4. For n < (d+2)/4 we should introduce appropriate counterterms. Of course, 
the one-loop effective potential may make perfectly good sense even when this Taylor series expansion is problematic. 
By this we simply mean that the one-loop contribution can be finite for values of the parameters and momentum that 
lie formally outside the radius of convergence of the above series: one would then obviously not proceed by expanding 
and integrating term-wise, as we have done here. 

Finally we should remind the reader that even if the generic reaction-diffusion-decay system is non-renormalizable 
for d > 4, this does not mean that such theories are completely useless. (Though it must be admitted that the 
number of physically relevant examples in four or more space dimensions is rather limited, [see for example variants 
on the idea of Kaluza-Klein theory], the non-renormalizability per se is not the issue.) All that non- renormalizability 
implies is that the theory must be viewed as an "effective field theory" that must include many more terms in the 
effective action than naively arise in the zero-loop approximation. These new terms carry with them additional 
(typically dimensionful) coupling constants, and these coupling constants govern the range of validity of the effective 
field theory. See jl3| for a modern discussion of effective field theories in the QFT context. In statistical mechanics 
language, the universality class of an effective field theory is much more complicated than would be naively deduced 
from the zero-loop approximation. 

We next plunge into a discussion of the above as a function of the spatial dimension, d. 



1 Since this point may cause some confusion to the reader, we briefly pause to belabour it: UV renormalizability is not the 
same as UV fimteness and we do not claim that the theory has to make sense at arbitrarily short distances. (If nothing else, 
in real chemical kinetics the interatomic spacing will provide a natural UV cutoff.) What UV renormalizability does is to 
sharply limit the number of relevant operators, so that the low energy theory (long distances, large times) is guaranteed to be 
relatively simple. Lack of UV renormalizability is not fatal for a theory, but does make life considerably more complicated. 
This phenomenon is known in the quantum field theory literature as "decoupling"; see e.g., Weinberg II ]13[ . 

2 See section (8.380.3). 
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A. Reaction- (diffusion)-decay: d — 



In d = space dimensions, there is of course no diffusion, and the system reduces to a simple Langevin reaction-decay 
system. We have 



df 



<t> = P{(j>) + V 



Specific examples of this behaviour include the noisy logistic equation (with additive noise) 



rd> 1 



d 

and the noisy Lotka-Voltera equation (used as a model for predator-prey interactions) 



(30) 



(31) 



—N = N(a + bP) + m , 

^-P = P(cN -d) + VP . 
at 



(32) 
(33) 



More generally, one can easily construct noisy versions of standard toy models as the Michaelis-Menten model for 
enzymatic autocatalytic reactions, the Goodwin switch (a model of feedback control), the Brusselator, the Fitzhugh- 
Naguno model of nerve potentials, or the Field-Noyes model for oscillating reactions jD. 

In d = there are tremendous simplifications in the general formalism. For time translation-invariant Gaussian 
noise it follows that 



V[0;0 o ;d = O] = ±P 2 (0) + ±.4 
Equivalently, 

V{fcfo;d = Q] = ±P 2 ( ( f>) + iA 



W) 



111 



1 



g 2 (m)P^)P"^) 
w 2 + [P'(0)]2 



4>o) + o(A 2 ) 



Alo 



(27T) 



In 



, 2 + [P'(«] 2 +j 2 ( W )PMF"W 



u 2 + [P'($)Y 



0(A 2 ) 



(34) 



(35) 



Just as in the case of field theory (d > 1) this is as far as we can go without making any further assumptions about 
the noise. For instance, temporally white noise implies 52 (w) ~> 1- Integrating over frequencies and using the integral 
( pp| ) supplemented by (^2|) and (p3|), yields 



V[4>\ ^0; d = 0] = iP 2 (0) + \A {R ev /[P'(0)]2 + p(0)P"(0) - v^PW } - (0 - <h) + 0(A 2 ) 
which can be re-written as 

V[4>; 4>o;d = 0] - iP 2 (0) + \A {ReyHPW 7 - \P'i,4>)\) ~ (0 ^ 0o) + 0(^l 2 ) . 



(36) 



(37) 



Note that from a field theory point of view a SDE in + 1 dimensions is almost quantum mechanics, and one-loop 
physics is almost semi-classical JWKB physics. To see what we mean by this, consider the quantum mechanics of a 
system governed by the Lagrangian 



('^) ~V(x) 



(38) 



The techniques more usually applied to quantum field theory can also be applied to quantum mechanics to obtain a 
quantum mechanical effective potential 



3M 



[x;x ] = !V(x) + \h Re 



V"{x) 



(X -> .T ) + 0(h 2 ) 



(39) 



This effective potential has the standard interpretation of being the minimum expectation value of the Hamiltonian 
operator, when extremized over stationary states satisfying (x) = x. In this case there is a second possible interpreta- 
tion in terms of the zero-point energy associated with the natural oscillation frequency f2 = y/V" /m, and it is in this 
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sense that one-loop quantum mechanics is equivalent to semi-classical quantum mechanics. The reaction- (diffusion- 
decay system is formally very similar to one-loop quantum mechanics with the replacement V(x) — > ^P(cf>) 2 ; the only 
difference arising from the manner in which the Jacobian is treated (the |P'(</>)| term). In quantum mechanics two 
classically degenerate minima often have their degeneracy broken via semi-classical effects: from the discussion above, 
analogous phenomenon is seen to occur in stochastic mechanics. 

(Note, however, that this is an analogy, not an identity. The SDE always gives rise to Wiener functional integrals, 
analogous to Euclidean quantum mechanics, instead of the Feynman functional integrals of quantum mechanics, and 
the SDE never exhibits the interference phenomena and complex amplitudes that are so characteristic of quantum 
mechanics.) 

B. Reaction-diffusion-decay: d = 1 

Spatial structures and patterns in one dimension, such as for example the prenatal tail markings in Genetta genetta 
(common genet) can be successfully modeled by one-dimensional reaction-diffusion equations. The reader is referred 
to the book by Murray for more concrete examples 0. 

For d = 1 the relevant integral, though finite, is not analytically tractable. We are interested in evaluating 

V[0;0 o ;d= 1] = iP 2 {0) + ^A J dk {lie^/[iyk 2 - P'(0)] 2 + P(0)P"(</>) - ^[uk 2 - P'(0)] 2 } 

-(<^0o) + O(-4 2 ) . (40) 
By expanding in a power series, rescaling, integrating, and summing we obtain the following result 

,, u , , , pa f , 1 , h^g \- (h\ F(l/2) T(2n - 3/2) ( P{cp) P"{cp) 

V[4>, 0o, d = 1] = \P (0) + -A vl/2 ^ y 2T(2n-l) { [P'm 2 

-(ct>^M + 0(A 2 ) . (41) 

This particular expansion only makes sense for P'(<fi) < 0, which is the region of field configuration space in which 
small perturbations die away in the absence of noise rl Furthermore, the radius of convergence of the resulting sum is 
equal to one, so that this expression is limited to the region PP" < (P 1 ) 2 - These limitations are not fundamental, but 
are artifacts of the expansion and integration procedure (the individual limits of first expanding and then integrating 
term by term do not commute) which must be taken into consideration. 

It is easy to see that the integral of equation fl4C|) is real, convergent, and well behaved for P'(4>) > and that 
P(4>) P"{(j>) > [P'(<fi)} 2 - Indeed, there is an exact (if rather formal) representation of this integral in terms of 
incomplete elliptic integrals of the first and second kinds. The integral can also be evaluated in terms of a 3F2 
generalized hypergeometric function and/or an assortment of complete Elliptic integrals Though exact, these 
formulations are too cumbersome to be useful and we can better understand the general situation by rescaling the 
original integral, using k — y/\P'(<fi)\/v x, to give 



V[0; d = 1] = |P 2 (0) + ±A |P 'ff 2 J™ dx jRe^il) 2 



PWP'W , 2±1| 



TO)] 2 j 

-(t^M + OiA 2 ) , (42) 
where ± = —sign[P' (</>)]. Thus without any detailed calculations we know that the form of the effective potential is 

-pf» p"{<t>r 



1 \P'((h)\ 3 / 2 
V[0;0 o ;d = 1] = hP 2 i<P) + tt A 12 F ± 



[p'm 2 



0o) + 0{A d ) , (43) 



with F± [z] the function 



3 Inserting <j>(t) = 4>o +S<j)(t) into the noiseless version of (^|), and expanding to first order in the time-dependent perturbation 
54>{t) yields: d[6<f>(t)]/dt = P'(4>a) S<j>(t) + 0[{8<f>) 2 \. So for P'{4>o) < perturbations are damped in the absence of noise. 

4 For further details see fiM, section (9.14), page 1045, and/or section (8.1). 
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F±[z] = J dx{Re v / (x 2 ±l) 2 +z - |.x 2 ±l|} , 



(44) 



such that F±[z = 0] = 0. It is that corresponds to the case P'{4>) < discussed above [see equation (|d|)]. The 

case F-[z] is trickier as there is no simple Taylor series expansion around z = (at least not in integer powers of z). 
The overall shape of these functions can be seen in Figures (pi) and (0) . 




FIG. 1. Plot of F±(z) for positive values of z, from zero to 100. 



F+ (z) and F- (z) 



or 



F+(z) 




F-(z) -1 -0.8 -0.6 -0.4 -0.2 

Z 



FIG. 2. Plot of F±(z) for negative values of z, from —1.00 to zero. 

Near a stable fixed point one has P'{4>) < 0, as mentioned above. We can use either equation (|4l]) or expand 
equation ( p4| ) for F+(2) to deduce that 

v 11 iP 2au 1 ^ T(j) 2 P^)P"(<p) ( P 2 (P"f \ , or , 2 , ,,,, 

This is the universal behaviour near a stable fixed point of the PDE. 

On the other hand, at a point of neutral stability for the original noise-free PDE we have P'(<fi) ~~ > 0> corresponding 
to z — > oo. Fortunately, the large z behaviour can be analyzed analytically. Some tedious integral analysis leads to 
the expressions 
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.l_„z 3 / 4 [l + 0(l/z)} 



F±(z) - ^ +oc 

-> 1.23605... z 3/4 (l + 0(l/z)) 



i r(|) 2 

6 



MS 



Z 3/4 (1 + 0(1/2)) 
Z — » +OO. 



P ± (-z) = z 3 / 4 [1 + 0(l/z)] - z 3 / 4 (1 + 0(l/z)) 

6 V27T 

-► -0.87401... z 3/4 (l + 0(l/z)) as z -► +00. 
Thus near a point of neutral stability we find universal behaviour described by 



(46) 
(47) 

(48) 
(49) 



V n outral[^; 00? d = 1] = ±P 2 (</>) + —-4 A± c 

Z7T 



[P(0)P"(^)] 3 / 4 

1/1/2 



O ) + O UP'(<^) 2 {/P(0)P"(0) + 0(A 2 ) . (50) 



Finally, near an unstable fixed point P'{4>) > 0. Then we must use equation ( |44| ) for P_(z). For z > a (not 
particularly obvious) series of manipulations leads to 



dF_(z) 
dz 



d.r 



l) 2 + z 



= iRe 



K 



1 



(51) 



Here K(fc) denotes a complete elliptic integral of the first kind. Unfortunately, although well-adapted to numerical 
work, this expression is not particularly illuminating from an analytic perspective and we do not pursue it any further. 

To summarize: we see that in d = 1 space dimensions a lot can be said about the structure of the effective 
potential, without ever having to specify the precise nature of the driving term. We trust that the general outline of 
the calculation is clear, and that applying the method to specific examples will now be straightforward. 



C. Reaction-diffusion-decay: d = 2 

In two spatial dimensions, reaction-diffusion equations have been extensively employed as models for the formation 
of patterns on animal coats (such as leopard spots), wing-marking patterns on butterflies, or fingerprint development 
in humans, etc. [SJ. 

For d = 2 the integral needed for the calculation of the effective potential is analytically tractable. After a change 
of variable, x = vk 2 , and the introduction of an ultraviolet cutoff, A, the desired integral becomes 

/ dx^(x - P') 2 + PP" = ±py (P') 2 + PP" - IPP" \n[^(P') 2 + PP" - P'] 
Jo 



+ i(w\ 2 - P')^/{vA 2 - P') 2 + PP" 
+ IPP" ]n[y/(vA 2 - P') 2 + PP" + vh 2 



(52) 



This integral is tabulated in reference ]lq| , equations (2.261) and (2.262.1). For the time being we assume that all 
arguments of both the square roots and the logarithms are positive, and leave the more technical details for later. 
The integral that appears in the one- loop effective potential, equation (|24|) , yields 



d.r 



VOr - P') 2 + PP" - V( x - P') 2 = h p W( p ') 2 + pp " ~ i( p ') 2 ~ i pp " HV( p ') 2 + pp " ~ p '\ 

+i(^A 2 -P') 2 



PP" 



{vk 2 - P') 



- 1 



\PP" ln(2i/A 2 ) 



\PP" In 



v /( i /A 2 -P') 2 + PP" + f A 2 - P' 



Taking the A — > oo limit, we get 



dx 



o 



y/{x - P') 2 + PP" - y/{x- P') 2 



IP' 



a/(P') 2 + PP" - P' 



2vA 2 



hPP" In 



pp" - p' 



(53) 



(54) 



- 1 PP" ln(2r/A 2 ) + i PP" + 0[l/(vA 2 )} 
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This explicitly verifies the presence of the logarithmic term expected from naive power counting |l9| . This logarithm 
is the only divergent contribution, and since it is proportional to P(</>)P"(0), the (one-loop) regularization may be 
performed by introducing the renormalization scale fx and making the following split into renormalized parameters 
and counterterms 

-Prenormalized (0) + A K P re normaIized 

(0) ln(A 2 /M 2 ) , (55) 

where if is a calculable numerical constant whose precise value is not important for the present discussion. It is useful 
to make an additional finite renormalization in order to eliminate the \PP" term from equation (|54]). (It is important 
to realise that these finite renormalizations do not affect the ground state structure. They are simply equivalent to a 
convenient choice of renormalization scale fj,.) 

After carrying out these steps, the one-loop effective potential becomes 



V[0;0 o ;d = 2] = iP 2 (0)+ - 4 



P' ( v/(P') 2 + PP" - P' ) - PP" In 



v/(P') 2 + PP" - P' 



-{ct>^<t> ) + O{A z ). (56) 

Like the action, the one-loop effective potential cannot explicitly depend on the renormalization scale \x. In fact, the 
renormalization group equation tells us that 

M-^V[0; O ; d = 2] = =>• = ~^P"{4>) + 0(A 2 ) , (57) 

which when combined with wavefunction renormalization and the general theorem of algebra will give the renormal- 
ization group equations for the couplings in V[4>; <po; d = 2] [pof . We note that this equation is similar in form to the 
one found for the P{4>)2 QFT in two spacetime dimensions, see p0[ . 

It is clear at this stage that the inclusion of a bare tadpole term (Po)bare is essential. If there was not a tadpole, 
then P(<f>) would start off as P\<\> + P2<f> 2 + Pa<f> 3 + • • *, so that the lowest order term in the zero-loop potential 
would be P 2 (f> 2 . On the other hand, as is explicitly seen in equation (|54|), the divergent terms are proportional to 

PP" = (Pi4> + P 2 4> 2 + P 3 </> 3 H )(2P 2 + 6P 3 </> H ) -> 2P 1 P 2 <j> H , and there is a divergent term proportional to cf> 

which is not present in the tree- level effective potential. Thus, in order to render the theory one- loop renormalizable 
a tadpole term must be included in the tree-level potential. 

Notice that if P(4>) is odd [P(— 4>) = —P((f>)], then the reaction-diffusion-decay is symmetric under the following Z 2 
symmetry transformation 

<p — ► —ip and T] — * — rj . (58) 

This prevents the generation of any of the even power monomial contributions to the polynomial P(4>), including the 
tadpole, and one never needs to introduce the tadpole at the tree-level. Notice however that this symmetry is not 
relevant to a realistic reaction-diffusion-decay system, since it excludes any two-body reactions (in fact all 2n-body 
reactions) . 

Finally, we mention what happens when one has to deal with one of the branch cuts that we have temporarily 
suppressed for simplicity of presentation. Although the intermediate stages of the calculation are algebraically involved 
(and messy), the ultimate answer is simple: take the real part of the expressions above. 



D. Reaction-diffusion-decay: d = 3 



An example in which prediction and observation of pattern formation in three dimensions has created much interest 
is in the field of nonlinear optical systems. Recently, three-dimensional reaction-diffusion equations (of the Swift- 
Hohenberg type) have been derived for degenerate optical parametric oscillators in which three-dimensional Turing 
structures and spatial solitons have been predicted to exist as stable structures |^l|. In addition, morphogenesis and 
structural development in embryos are examples of intrinsically three-dimensional phenomena. 

The case d — 3 is a straightforward generalization of the d = 1 result. We are interested in evaluating 

Vfo; d = 3] = iP 2 (0) + j^A J™ dk 4irk 2 {liey/[isk* - P'(0)] 2 + P(0)P"(0) - VH 2 - P' (</»)] 2 } 

-{cfj^^ + OiA 2 ) . (59) 
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The integral is no longer finite and a single renormalization (without running logarithms) must be performed to absorb 
the infinity in the renormalized parameters. One must be careful to keep track of all the numerical coefficients and 
the final result is 



V[0;^>;d = 3] = iP 2 (4>) 



1 



[-P'{<t>)f 2 



(2nr 
) + 0(A 2 ) 



,3/2 



E 



L\r(3/2)r(2n-5/2) ( P{4>) P"(<f>) 



2r(2n-l) 

n— 1 x / v 7 



(60) 



This expansion again only makes sense for P'(4>) < (the region of field space in which small perturbations die away 
in the absence of noise). Furthermore, the radius of convergence of the resulting sum is one, so that this expression 
is limited to the region PP" < (P') 2 . These limitations are again not fundamental, but are merely a reflection of our 
choice of Taylor expansion in terms of the field variables. In a manner similar to what was done in d = 1 we can also 
write the above equation (60) as 



V[0;0 o ;d = 3] = iP 2 (0) + 



with P^ _3 [,z] now being the function 



(2n) 



JPW/ 2 

,.3/2 *± 



P(cj>) P"(0) 



TOP 



<h) + o{A 2 



F: 



rf=3r 



dx x 



He^{x 2 ± l) 2 + z 



\x'±\\~ h - 2 



(61) 



(62) 



again such that F±[z = 0] = 0. It is F+[z] that corresponds to the case P'(4>) < discussed above [see equation (p0[)] . 
The case P_ [z] is again trickier as there is no simple Taylor series expansion around z = 0. Note that the last term 
z/x 2 is the counterterm introduced to guarantee UV fmiteness, and that the integral is also IR finite. These functions 
can be analyzed in a manner analogous to the discussion for d = 1 but for the sake of brevity we do not repeat details 
which are left to the industrious reader. 



E. Special cases: monomial interactions 



There are nice simplifications for monomial interactions, where P{<j>) = £4> n . (Because of the generic presence of 
the tadpole term this monomial behaviour should always be imposed on the renormalized interactions, not the bare 
ones.) For monomial interactions the combination P(0)P"(</>)/[P'(0)] 2 reduces to the constant (n — 1)/"-, and the 
one-loop effective potential (for d = 1, and suppressing <f>Q for convenience) becomes 



V[0;d = !] = *£ 



2 ,2n 



1 J-£0 W ~ 

2n 1/I/2 



113/2 



-K ± [l;n]+0(A 2 



(63) 



with i£±[l;n] a calculable (n-dependent) dimensionless constant. If the differential equation (without noise) is 
assumed to be stable against small perturbations, the coupling £ must be negative, and the exponent n must be an 
odd integer n = 2m + 1. In this case P'{<j>) = (2m + l)£</> 2m < O ^ < guarantees that linear time-dependent 
perturbations will decay in time. The effective potential becomes 



vr 



1 \f\ 3 / 2 \(h\ 3m 

:d = l] = ^ 2 ^ 4m + 2 + —A ™ K+[l; 2m + 1] + 0(A 2 ) 



(64) 



From the previous equation, one can see that noise induced corrections to the effective potential dominate for small 
fields in one space dimension. 

For two space dimensions the restriction to monomial interactions implies 



Vty, d = 2] = + -^-£ 2 ^2n-2 I n f ^ 2n 2 _ - _ n )_ n t n _ x \ ln 



^ n_1 (V2n 2 -n-n) 



By making a finite renormalization fx — * (a, we can simplify the previous equation to obtain 



V[cj>\<t> ;d = 2] - ^ 2 4> 2n 



A 



i-2 ±2n-2 



l)ln 



n-1 



9 



0{A 2 



0{A 2 ) . (65) 



(66) 
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The coefficient of the logarithm is negative, implying a breakdown of perturbation theory for large fields. More 
importantly, since the coefficients of P(<j>) run at one-loop, and a monomial is not a fixed point of the renormalization 
group equations, if we tune the interaction to be monomial at some fixed scale p, then (in d = 2) the interaction will 
not remain monomial if the scale is changed. 

For three space dimensions the situation is similar to d = 1. For a monomial interaction the one- loop effective 
potential is given by 

1 I _ C ^n-l 15/2 

V[0; O ; d - 3] = if 4> 2n + j2nV A - J3/2 K ± & n \ + °( A ) • ( 67 ) 
For a system stable in the absence of noise (£ < and n = 2m + 1), we have 

V[</>; O ; d = 3] = K 4m + 2 + 7^2 -A - ^ K + [3- 2m +1} + 0{A 2 ) . (68) 
In this case the noise induced effects (one- loop effects) become important for strong fields. 



F. Special cases: higher dimensions 

We have seen that arbitrary polynomials in the reaction-diffusion-decay system subject to white noise are one-loop 
ultraviolet renormalizable in 0, 1, 2, and 3 space dimensions (in fact, finite for and 1 space dimensions). We can 
extend the range of dimensionalities in which these systems are one-loop renormalizable at the cost of restricting the 
form of the interaction. 

For instance, in four space dimensions we have already seen that there is a divergence proportional to 
P(4>)P' '(<f)) P" ' ((f)). For general P{(f>) this cannot be renormalized, but if P(4>) is a polynomial of degree three or less, 
then P' \4>)P" '(</>) is also a polynomial of degree three or less. For this restricted class of interactions, the divergence 
can be absorbed into bare potential even in four space dimensions. 

In six space dimensions there are two new divergent terms. They are proportional to P(<f>)[P' (</>)] 2 P" (4>) and 
[P((f))P h '(4>)] 2 , respectively. If P(4>) is a polynomial of degree two or less then [P 1 ((f))] 2 P"(4>) is also a polynomial of 
degree two or less and the renormalization program can be carried out. In this case P"((f>) is either a constant or 
zero, so that the second type of divergence is no further obstruction. 

Finally, in eight space dimensions there is only one new divergent term. It is proportional to P(<fi)[P' ((f))] 3 P" ((f)) 
and so the theory is one-loop renormalizable only for linear interactions (i.e., for free fields where the theory is not 
only renormalizable but is actually finite.) 



G. More general noise 

It is clear from the above arguments that the ultraviolet renormalizability of the reaction-diffusion-decay system 
depends critically on the large momentum behaviour of the noise two-point function. In some problems, colored noises 
maybe of interest, and we just give a necessarily brief discussion of the modifications needed in our analysis. Let us 
suppose that the noise is more general than (space-time) white noise. For instance, let us assume the noise is still 
temporally white, but spatially power-law distributed in the ultraviolet region with 

92$) = h(\k\) « (k/k )- 9 6(A - k) , (69) 

where the positive exponent 9 characterizes the strength of the ultraviolet-singular noise. In order to obtain the 
divergence structure of the one-loop effective potential, we must make use of equation (|24|). It is then easy to 
see that the first two terms in the expansion for the effective potential have ultraviolet behaviour proportional to 
P((j))P"((f>) A d - 2 - 6 , and P((j>) P'(4>) P"(4>) A d - 4 ~ e , respectively. Since the bare potential is P 2 (</>), to have any hope 
of absorbing the infinities into the bare action we must have d < 4 + 0. The one-loop effective potential for the 
reaction-diffusion-decay system is then one-loop ultraviolet renormalizable for d < 4 + 9 spatial dimensions. 



If the noise is not temporally white, (but still Gaussian) , one must revert to equation ( 18 ) and perform a case-by-case 
study. 
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III. NOISE AND STABILITY: AN APPLICATION 



In this section we illustrate in broad strokes how the one-loop effective potential can be used to investigate the onset 
of instabilities and pattern formation in physical, biological, and chemical systems modelled by reaction-diffusion 
equations 

By way of a concrete example which will serve as a simple template for the one-loop equation, we consider the 
following model: 

^-vV 2 <j> = P(<j>)+r,(x,t), (70) 
at 

where P(<f>) = a<t> 2 + bcf> + c is the kinetic reaction polynomial parameterized by three real constants a, b, c, and the 
diffusion constant v is a positive real number. We take the noise to be Gaussian and white. We first briefly run 
through the standard steps needed to perform a linear stability analysis of the noiseless or zero- loop version^ of (|7(]) . 
This will serve as a point of reference when we come to discuss the linear stability analysis to be performed on the 
effective one- loop version of (f70|). 



A. Zero-noise analysis 

As is well known, the study of the onset of symmetry breaking instabilities, be they Hopf bifurcations or Turing 
instabilities, starts by classifying all the static and spatially homogeneous solutions <fio of the reaction-diffusion equation 
at hand. These constant field configurations represent the maximally symmetric states of the system, which could be 
stable or unstable with respect to time and/or space dependent disturbances. For our toy model, these states satisfy 

P(<t> ) = , (71) 

which is the result of evaluating (|7^) for constant fields in the absence of noise (which is completely equivalent to the 
zero- loop limit of this equation) . The solution is immediate, namely we have that 



, -b ± V6 2 - 4ac 
4 = Ya ■ (72) 

It should be pointed out that if the field <fi represents a chemical concentration, then (f> and 4>q must be non-negative. 
There are many choices of the control parameters a,b,c for which this condition is met. With knowledge of the 
constant states, the next step is to expand about them so that we can study both the temporal and spatial evolution 
of disturbances with respect to these symmetric states. Setting ip^ — <ft — </>q , leads to 

which is the exact zero-noise equation for the perturbations. In arriving at this expression, we have used ( |7l| ) and 
P'((/)q) = ±\/6 2 — 4ac, P"((f>Q) = 2a. To study the onset of linear instabilities, we need only focus on the linear part 
of this equation. This is most conveniently carried out in momentum or mode space, for which we introduce the 
Fourier transform (in any number d of spatial dimensions) of the fluctuation ip(x, t) 



vV 2 (p± = ±cp± V& 2 - 4ac+ a(y ± ) 2 , (73) 



/ Td^'^W- ( 74 ) 



d d q 

In terms of the mode functions <p$(t), the (exact) zero- noise equation for the perturbations takes the form 

d d k _ _ 



dt 



± \/b 2 — 4ac iff + a J 



(2tt) 

= X(q 2 ) (f>q-+ non — linearities , (75) 
which identifies the momentum-dependent eigenvalue X(q 2 ) and the non- linear mode-mode coupling terms. 



5 The noiseless limit of general reaction-diffusion equations need not necessarily coincide with the zero-loop limit, though these 
limits are in fact identical for our model equation; see the further comments to this effect below. For general details regarding 
the distinction between no-noise and zero-loops, see reference J7|. 
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1. Onset of Hopf bifurcations 



By definition, Hopf bifurcations are spatially homogeneous but time dependent linear instabilities jjj. Spatial 
homogeneity corresponds to the zero- momentum mode (q = 0), so the onset of this class of instability is revealed 
by studying the properties of the zero- momentum eigenvalue: X(q 2 = 0). A Hopf instability occurs when the real 
part ReA(O) = 0. It is useful to consider the mode eigenvalues as functions of the control parameters (a, b, c) in the 
reaction kinetics polynomial. If we define OHopf by b^ op{ = 4ac, then it is easy to see that for b < &Hopf , the eigenvalue 

A(0) is pure imaginary and the linearized perturbations oscillate about the constant and static field configurations <pQ 
and (j>Q with the same frequency u> = \b 2 — 4ac\^ . On the other hand, if b > &Hopf, then the linearized perturbations 
about (f>Q decay exponentially, and they grow exponentially about 0j. The instability therefore sets in at b — fefjopf, 
and only the <pQ field configuration is unstable. 



2. Onset of Turing instabilities 

The onset of spatial structure formation (Turing instabilities) occurs whenever the eigenvalue satisfies the condition 
Re[A(q 2 )] > 0, for some non-zero mode or modes 0]. The corresponding mode (or modes) set the length scale L (or 
scales) that characterize the spatial structures: L ~ 1/q. Note that for b < &Hopf, the real part of the eigenvalue for all 
modes is negative, so no spatial structure can form, and this holds for both the initial configurations (jy^ . Furthermore 
for <pQ , b > &Hopf, and for any q, the eigenvalue X(q 2 ) is a negative real number. Therefore, no spatial patterns can 
develop as linear perturbations of this field configuration. 

On the other hand, consider the case 4>o an< i b > &Hopf • The eigenvalues for all modes are real, and of these, there is 
a finite band of momentum scales for which the eigenvalue is strictly positive: namely, for < q 2 < q\ = y/b 2 — 4ac/v. 
So, one expects onset of spatial structures to form with length scales corresponding to this momentum band. 

This completes the purely linear stability analysis of the toy reaction-diffusion equation ( |70| ) in the absence of noise. 



B. One-loop analysis 

We now demonstrate, (by making use of the effective potential calculated in previous sections), how the inclusion 
of noise at one-loop impacts on the linear stability analysis illustrated above for our simple reaction-diffusion model. 
The idea is to repeat the above steps, but working now with the effective one-loop reaction-diffusion equation. To 
obtain the equations of motion in the presence of noise, we refer to the discussion in Appendix B of |t]]. The complete 
effective action for reaction-diffusion systems in the presence of (arbitrary) Gaussian noise is given by 

<M = \j Jd d xdt d d ydt' (D cS cp - P cS [4>])g^{D cS cp - P eff [<f>]) , (76) 

where D e g and P e ff are an effective differential operator, and effective kinetic reaction function, (not necessarily a 
polynomial!), respectively. The full (to all loops) equation of motion follows from the condition jffl: 

^ = - ^ = 0, (77) 
00 dtp 

where the second expression holds for the effective potential and yields the dynamical equation satisfied by static 
and homogeneous field configurations. As a quick check of the formalism, consider the zero-loop effective potential: 
V[0; 4> ] = iP 2 [0]- From (Q) we must have P(0)P'(0) = 0, so that either P(<f>) = and/or P'(0) = 0. In fact, we 
know from ( fn] ) that it must be P(</>) that vanishes to correctly yield the solutions of (JT0|) in the noise-free static and 
homogeneous field limit. The "spurious" factor P 1 ((f) is a consequence of the quadratic nature of the effective action 
([76]). It is easy to verify that P'(4>) ^ evaluated at the zeroes of P((fh This is important in that we do not generate 
more solutions than those corresponding to the zero-noise equation (|70|): zero- loops should correspond to zero- noise, 
and we find that this does in fact hold for our model. In general, spurious solutions will be absent whenever a certain 
differential operator [D — P'[(j)}) is non-singular (invertible) . Even in those situations where it can become singular, 
a simple limiting procedure can be invoked to eliminate the spurious solutions Q . 

Under these conditions, i.e., making use of the invertibility of (D— P'[(f>]), the one-loop equation of motion associated 
with (0) is 
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^- ! /VV + O(^lV 2 0,^,---) = (2V[0;0o])' = [P 2 (cf>) + 2AX(cb) + 0(A 2 )]> , (78) 
where the one-loop terms in the effective potential denoted above by X(<p) have been computed for general kinetic 



reaction functions and may be written down by inspection from (|37|), (43), (pq), and ( pCf ) above. It should be noted 
that the one-loop contribution to the effective potential vanishes whenever P((jr) = so that without loss of generality 
we can write the O(A) terms as 

X{4>) = P{4>) h{4>) , (79) 

with h(<f>) remaining finite (and generally non-zero) as one approaches solutions of the zero-noise equations of motion. 
Furthermore, at one-loop we know from (rather general considerations detailed in fll9| ) that the one- loop effective 
action for any reaction-diffusion equation driven by white noise has no wavefunction renormalization in fewer than six 
spatial dimensions. This means that in dimension six or less, the differential operator D per se remains unchanged 
at one-loop. There can, however, be finite renormalizations that induce finite and calculable one-loop corrections 
that involve combinations of fields and derivatives; we summarize the first few possible structures belonging to this 
class above. It is clear from (|7^) that as the noise amplitude A is taken to zero, we recover the zero-loop (noiseless) 
equation (|70|) of motion. 

In view of the above, the one-loop equations of motion read 



g - vv 2 4> + 0{AV 2 4>, Ad 2 4>, ...) = P{4>) 



P(cj)) + Ah{(t>) + 0{A 2 ) . (80) 



We are now ready to proceed with the stability analysis for the one- loop reaction-diffusion equation (|78|) . Just as for 
the noiseless, zero-loop case, one begins by solving for all the possible static and spatially homogeneous configurations. 
These will be the solutions of 

P(M+Ah(4, a ) + 0(A 2 ) = 0, (81) 

no matter how complicated the derivative structure on the left-hand side of the one-loop equations of motion may be. 
We denote by 4>q the one-loop constant field configurations to distinguish them from their tree-level counterparts. In 
fact, it is easy to demonstrate that the zero-loop solutions (f>o are in general not solutions of this one-loop equation 
©. Thus, we can conclude that typically (j>Q ^ fa. Physically, this reflects the fact that the presence of noise 
has altered the symmetric states of the system. Suppose we have chosen a particular spatial dimension and have 
catalogued these new symmetric states. [To do so in actual practice requires selecting one of (|37|), ([f3"|), (|56|), and 
(|60|), and solving the resulting (algebraic-transcendental) equation implied by (^). However, bear in mind that the 
point we wish to make can be achieved without doing so explicitly] The next step involves expanding the one-loop 
equation in (linear) perturbations about these one-loop states. We define ip = <f> — (f>o, and write 

% - V V 2 + 0(AV 2 & AdU, ...) = P(ip + 0o) + A h(<p + 4> ) + 0{A 2 ) . (82) 
at 



We Taylor expand the right-hand side of ( J82| ) making use of ( J8l| ) . We then obtain the linearized one-loop equation of 
motion 

- vV 2 (p + 0(AV 2 tp, Ad 2 ip, ...) = \p'(4>o) + A h'Uo)] <t> + 0(A 2 ) + non - linearities . (83) 
ot L J 

After Fourier transformation the evolution of the mode functions is governed by 

^ = A( q 2 ) ^ + 0{A 2 ) + nonlinearities - 0{Ad 2 t ~^^ . . .) , (84) 
where the one-loop mode-dependent eigenvalues are 

\{q 2 ) = -vq 2 + P'($ ) + Ati(4> ) - 0(Aq 2 , . . .) . (85) 

The function h is the one-loop contribution to the effective potential (divided by P) and has been explicitly calculated 
(for arbitrary P) in various spatial dimensions. Since it is a function of the reaction polynomial, it is also a function of 
the control parameters (a, b, c) appearing in P. If there arc finite renormalizations leading to new derivative structures, 
then A will also depend on their respective numerical coefficients. For the toy model considered here, this means that 



1G 



the one-loop mode eigenvalues A are also functions of the same parameters (a, b, c) that appeared in the zero-loop 
eigenvalues in (|75|), as well as of whatever finite renormalizations are present at one- loop. 

We summarize our main point: when we investigate the conditions for the onset of Hopf and Turing instabilities 
at one-loop, Re[A(0)] = and Re[A(q 2 )] > 0, q 2 7^ 0, we will obtain new conditions on the parameters that govern 
the onset of whatever instabilities are present, with respect to the 0(A) noise- altered symmetric states <f>Q. We see 
that the effect of the noise is to shift the symmetric states of the system, as well as to change the nature of the linear 
instabilities that may be induced by perturbations around these new states. 



IV. DISCUSSION 



In this paper we have made use of the formalism developed in jTj and have applied it to the class of stochastic partial 
differential equations subject to diffusion, decay, and polynomial reaction kinetics. These equations are, for example, 
extensively considered in the mathematical modeling of chemical reactions, and biological pattern formation 

We have restricted that framework to the case of white noise and have calculated the one- loop corrections (i.e., we 
have taken into account the second-order fluctuations about the static and homogeneous solutions of these equations) 
to the effective potential in various spatial dimensions and for a general polynomial reaction kinetics term P((j>). The 
effective potential, which provides information about the possible ground states of the system (which may or may not 
be stable to small perturbations, see below), is calculated by functional integral methods and we find that it is one-loop 
finite for zero and one space dimensions, and one-loop renormalizable in two and three space dimensions. By finite we 
mean there are no short-distance divergences, and by renormalizable, we mean that whatever ultraviolet divergences 
are present, they can be absorbed into the parameters appearing in the original stochastic partial differential equation. 
In particular, in two space dimensions, this renormalizability leads to a set of one-loop renormalization group equations 
( p7| ) that govern the renormalization-scale dependence of the parameters present in the original reaction polynomial. 

In understanding the onset of spatio-temporal pattern formation in systems out of equilibrium, it has proven 
extremely useful to begin the analysis by first solving for and then classifying all the static and spatially homogeneous 
states allowed by the time-dependent partial differential equations employed to model the system in question . In 
this way, one can straightforwardly decide whether the system will exhibit Hopf bifurcations and/or Turing instabilities 
and get a handle on the qualitative nature of the pattern expected to emerge. 

This can be followed up by an amplitude analysis of the fluctuations about these static and homogeneous states. 
The unstable modes are the ones that lead to non-trivial patterns. For out-of-equilibrium systems coupled to noisy 
environments (or with inherent internal noise) it is important to know how the stochastic sources can alter and shift 
these static and homogeneous states, since these affect the onset of the pattern-forming (linear) instabilities. It is 
seen that indeed noise will affect these patterns in a way which is computable. In fact, the calculations in this paper 
show how the effects of stochastic noise on these states of reaction-diffusion systems can be taken into account in an 
elegant and computationally direct way following the general formalism developed in H . 
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APPENDIX A: JACOBIAN FUNCTIONAL DETERMINANT IN ZERO SPACE DIMENSIONS 



We are interested in evaluating the following functional determinant in zero space dimensions 

dP" 

J = det (d t ) n -- 
d< 

If n > 1, from the general analysis given in 0, we see 



(Al) 
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J n = det(ft)" , 



(A2) 



whereas for n = 1, we get |?]] 



dP\ 



Ji = det(ft) exp ( -6(0) - -> dot (ft) exp 



1 dP\ 

2 d^J 



(A3) 



where we have adopted the prescription 6(0) = 1/2. 



APPENDIX B: JACOBIAN FUNCTIONAL DETERMINANT FOR THE 
REACTION-DIFFUSION-DECAY SYSTEM 



For the reaction-diffusion-decay system F[<fi(x, t)] = P(<ft) is an ultra-local function of </>, most commonly a polyno- 
mial in d>. Thus 



SF[</>(x,t)] 
84>(y,t') 



P'{4>(x))5(x-y)8{t-t') 



The relevant trace is given by 



Tr 



Gi- 



84>(y) 



6(0) / ditr spaco [P'{4>) S(x-y)} = 6(0) S d (0) / d d x dt P'{4>) . 



(Bl) 



(B2) 



In contrast to the Kardar-Parisi-Zhang equation jl7| , the Jacobian determinant for the reaction-diffusion-decay system 
is not a field independent constant, but (adopting the prescription 6(0) = 1/2) one has 



Jrdd = exp 



±6 d (0) ! d'.F.l/ />'(„) 



(B3) 



Thus for the reaction-diffusion-decay field theory the functional determinant at worst leads to extremely simple 
Faddeev-Popov ghosts. There are general arguments (see for example Zinn-Justin jl(| pp. 373, 307, or related 
comments in Itzykson-Zuber 15 p. 448) to the effect that terms proportional to 5 d (0) can always be safely discarded 



in dimensional regularization. We do not want to step into the middle of this contentious issue and merely note that we 
have found it more convenient to not adopt the formal result S d (0) — 0, and instead to explicitly carry the Jacobian 
along in the calculation. Keeping the Jacobian explicit is essential to showing one-loop finiteness in d — 1 space 
dimension, a result that would otherwise be disguised by unnecessarily discarding the precise counterterm needed to 
ensure the one-loop finiteness of the theory. 

This situation is in marked contrast to that for the KPZ system [|T7|j , wherein the Jacobian determinant is a field 
independent constant irrespective of how one wishes to handle the formal result S d (0) = 0. 



APPENDIX C: FEYNMAN RULES FOR THE REACTION-DIFFUSION-DECAY SYSTEM 



We have derived the general form of the Feynman rules in the direct formalism, applicable to arbitrary SPDEs, 
m @. There are a number of technical simplifications for the reaction-diffusion-decay system which make it worthwhile 
to present this particular case. From the reaction-diffusion-decay stochastic differential equation (ft) 



d_ 

dt 



v{V 2 - m 2 ) 



= p{<t>) + n 



we deduce the characteristic functional (partition function) Q 



Z[J] 



(t>4>) Vj J f 



exp 



ft^- i/(V 2 - m 2 )(j) - P((j)) 



x exp 



[J4>] 



(CI) 



d t (j) - y{V 2 - m 2 )(j) - P((j>) 



(C2) 
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from which we can immediately deduce the associated Feynman rules. Note we have opted to collect whatever terms 
linear in <j> there may be in the reaction polynomial P((f>) and place them on the left hand side of the equation. 
Compare to the equations (|^) and (|ll]) in the introduction: the decay rate is simply 7 = vm 2 . The functional 
determinant can cither be calculated from the preceding section, or it can be included as a Faddeev-Popov ghost 
term, with the same effect. There is only one field propagator and two vertices. The propagator is identical to the 
KPZ case, but the vertices have as high an order as determined by the highest power in P (see below). The following 
holds for translationally invariant noise. 

These are the basic structures (in the interacting classical action) from which one can derive the appropriate vertices 
by taking the corresponding functional derivatives: 



G ficld (fc, w) = - G ^ u) — , (C3) 



2 



5(k 1 + k 2 ) %i + wa) , (C4) 



Propagator 

;,u1 = 

lu 2 + v 2 {k 2 + TO 2 ) 

P(4>)-(j> vertex: 

[P(0)](fci,^i) 4>(k 2 ,uj 2 ) [-iui + v(k\ + to 2 )] 
(2^) rf+1 Wl ) 
P(4>)-P(<t>) vertex: 

1 [P(0)Pi,^i) [P{4>)]{k 2 ,u 2 ) g g 

- ■ =; o{k 1 + k 2 ) d(LOi+uu 2 ) . (05) 

2 (2^) d +iG r) (fc 1 ,c 1 ) 

For a polynomial P(<j>) of order n, there will be up to 2n propagators meeting at a vertex. Calculations with these 
Feynman rules involve fewer propagators and vertices than those of the Martin-Siggia-Rose formalism The 
basic trade-off is this: one reduces the number of propagators and vertices at the cost of making the vertices more 
complicated. As long as one is interested in one-loop physics this cost is not prohibitive , and the pay-off in terms 
of ease of calculation for the effective action and effective potential is considerable. 

APPENDIX D: AN INTEGRAL 

We want to demonstrate that 

r +oc (u?-z 2 ±it\ , 

J M dUjln [ co 2 + Y 2 )=^^ Z - y ) ■ (Dl) 

To see this note that 

+ °° /, .2 <72_|_„v\ r + °° / ' l LU 2 — Z 2 \\ f +Z f\rj 



^ u2 + y2 ) = J oo ^—J.^]^ to-J^ d.ln^^ T ^j ± 2^. (D2) 
The remaining integral can be found in |18[ equation (2.736.1) and (2.733.1), and we can finally write 



2 v2| \ r - » v.' \ -/ - 1 +°° 



) +^ln 


oj + Z 






uj — Z 



2Y" tan 



■'(f) 



-2ttY . (D3) 



APPENDIX E: THE ITO CALCULUS VERSUS THE STRATONOVICH CALCULUS 

In evaluating the Jacobian functional determinant one encounters a factor of 6(0), which is ill-defined and must be 
specified by some prescription. The prescription which is most useful in this context, and which we have adopted in the 
bulk of this paper, is the symmetric prescription wherein 0(0) equals i. This may be justified by a limiting procedure 
as described, for example, in the text by Zinn- Justin |l(| (Chapter 4, pp. 69-70.) The symmetric prescription is 
equivalent to adopting the Stratonovich calculus for stochastic equations. Choosing O(0) = is equivalent to the 
Ito calculus. The Ito calculus simplifies the Jacobian determinant (to unity) at the cost of destroying equivariance 
under field redefinitions (the Ito calculus explicitly breaks coordinate invariance in field space). See, for instance, 
Eyink [O, or Zinn- Justin |16|. In this Appendix we sketch the modifications required to implement the Ito calculus. 
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These changes are straightforward if at times tricky (the loss of reparameterization invariance under field redefinitions 
implies that all arguments involving a change of variables must be carefully re-assessed) . 

If we adopt the Ito calculus, then for any SPDE the (unrenormalized) expression for the one-loop effective potential 
simplifies to 



Vito[^;^o] = hF 2 m + $A I ,,,, In 



D\k,uj) - 




D{k,uj) - 



SF 



+ g 2 (k,Lu)F[. 



S 2 F 



(27r) d + ] 

- ^ cj> Q ) + 0(A 2 ) . (El) 
[Compare with equation (|j).] When specialized to the reaction-diffusion-decay system this further simplifies to 

uj 2 + \vk 2 - P'(cf>)} 2 +g 2 (k,w)P( ( f>)P"(4>) 



u 2 + [vk 2 - P'(0 O )] 2 + g 2 {k, uj)P(MP"(M 



0{A Z ) . (E2) 



[Compare with equation (^9|).] Further restricted to temporally white noise we obtain 



d d fc 
(2^p 



Re 



^fc 2 -P'(</>)] 2 +.9 2 (fc)P(0)P"(^) 



0o) + 0{A 2 ) 



(E3) 



[Compare with equation (£4|).] 

On the one hand, this looks like a tremendous simplification of equation On the other hand, the ultraviolet 

renormalizability properties are now considerably worse. For instance, by inspection of the above it is easy to see that 
one can no longer rely on a cutoff in the noise spectrum, g 2 (k), to keep the effective potential finite. Instead a cutoff 
in the momentum integral must be introduced by hand, complicating the process tremendously. Even if g 2 (k) — * 
for large momenta, for d > there exists an UV divergence proportional to A d [P'((f)) — P'((j>o)]. For generic P(<f>) 
this cannot be absorbed into a counterterm in the zero-loop effective action. That is, adopting the Ito calculus for 
the reaction-diffusion-decay system results in a theory that is one-loop non-renormalizable for d > 0, and so must be 
viewed as an "effective field theory" . 

Thus even for d = 1, where the Stratonovich calculus leads to a one-loop finite result, the Ito calculus leads to 
complicated expressions which obscure the underlying physics. It is for this reason, (the complications of dealing with 
non-renormalizable effective field theories), and the fact that the Ito calculus is not invariant under field redefinitions, 
that we have not further explored the Ito calculus in this paper. 

The one case where the Ito calculus gives a simpler answer than the Stratonovich calculus is for d = 0, corresponding 
to stochastic mechanics rather than stochastic field theory. In that case 



6 ;d = 0] = ±P 2 {0) + ±A ReV^ 2 W>)] 



<fo) + 0(A 2 



(E4) 



Of course, one should not be alarmed that the Ito calculus and the Stratonovich calculus give different intermediate 
results — they are different theories. Because the Ito calculus is not invariant under field redefinitions it is generally 
possible to find some choice of field variables (a reparametrization) that makes the two systems agree with each other, 
but that special set of field variables may not be the ones one naively started out with. 
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